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摘要 :【 目的 ] KEAS k Leguminivora glycinivorella ( Matsumura) 是 我 国 大 豆 Glycina тах 最 重要 的 
星 英 害 虫 。 本 研究 旨 在 探讨 中 国 大 豆 食 心 虫 不 同 地 理 种 群 的 遗传 分 化 情况 。【 方法】 通过 测定 大 
豆 食 心 虫 19 个 地 理 种 群 208 头 老 熟 幼 忠 的 线粒体 CO 开 基 因 序 列 (723 bp) ,利用 MEGA6.0 和 
DnaSP 5.0 等 软件 分 析 大 豆 食 心 虫 不 同 地 理 种 群 的 基因 序列 变异 和 遗传 分 化 情况 。 【结果 】 在 19 
个 大 豆 食 心中 地 理 种 群 中 共 获 得 208 条 CO 开 基 因 序 列 ,发 现 61 个 变异 位 点 ,包含 57 个 单 倍 型 。 
总 群体 单 倍 型 多 样 度 Hd =0.8522 ,各 地 理 种 群 单 倍 型 多 样 度 范围 在 0.1818 ~ 0.9546 间 ,总 群体 的 
ХАЖ Fst =0.4619, 总 遗传 分 化 系数 Gst =0.1186, 总 基因 流 Nm=0.29。 总 群体 Tajima”sD 检 
验 结果 不 显著 ,表明 大 豆 食心虫 在 较 近 的 历史 时 期 未 出 现 群体 扩张 现象 ,群体 大 小 稳定 。 单 倍 型 网 
络 图 、ML 单 倍 型 系统 树 及 采用 不 同 地 理 种 群 间 遗 传 距离 构建 的 UPGMA 系统 发 育 树 均 显 示 , 贵阳 
(GY) 和 都 安 (DA) 种 群 与 其 余地 理 种 群 分 化 明显 。AMOVA 分 子 变异 分 析 结 果 显 示 , 东北 组 .西北 
组 、 黄 淮 组 、 华 南 组 及 西南 组 种 群 组 间 所 占 总 变异 的 比例 大 于 组 内 。【 结论 】 大 豆 食 心 虫 不 同 地 理 
种 群 间 基 因 交 流 较 少 ,整体 的 遗传 多 样 性 和 遗传 分 化 程度 较 高 ,遗传 分 化 主要 来 自 地 理 种 群 组 间 ， 
且 地 理 距离 是 影响 不 同 地 理 种 群 间 遗 传 距离 的 重要 因素 。 
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Analysis of the genetic differentiation among geographic populations of 
Leguminivora glycinivorella ( Lepidoptera: Olethreutidae) in China based 


on mitochondrial СОП gene sequences 

ZHU Shi-Yu, XU Wei, САО Үп, CUI Juan, OU Shi-Qi, SHI $һи-беп * ( Innovation Center of Soybean 
Region Technology, College of Agriculture, Jilin Agricultural University, Changchun 130118, China) 
Abstract: [ Aim] Leguminivora glycinivorella ( Matsumura) is опе of the major pests of soybean Glycina 
max in China. This study aims to investigate the genetic differentiation among different geographic 
populations of L. glycinivorella in mainland China. [Methods] A 723 bp segment of the mitochondrial 
cytochrome с oxidase subunit П ( СОП) gene in 208 mature larvae from 19 geographic populations of L. 
glycinivorella in China was sequenced and analyzed. The genetic differentiation and sequence variation 
were analyzed using MEGA6.0 апа DnaSP 5. 0. | Results] Among the alignment of 208 СОП sequences 
from 19 geographic populations of L. glycinivorella, 61 variable sites and 57 haplotypes were identified. 
The total haplotype diversity index ( Hd) was 0. 8522, and the range of haplotype diversity of different 
populations was 0. 1818 — 0. 9546. The total fixed coefficient ( Fst), total genetic differentiation 
coefficient ( Gst) and total gene flow (Nm) меге 0.4619, 0.1186 and 0.29, respectively. No significant 
difference among different populations was found by Tajima’ s D test, suggesting that there might be no 


population expansion in recent history of L. glycinivorella and the population size is stable. Analysis of 
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haplotype networks, МІ, phylogenetic trees of haplotypes апа UPCMA phylogenetic trees based on the 


genetic distance of different geographic populations showed that Guiyang (СҮ) and Du’ ап (ПА) 


populations had obvious genetic differentiation with other populations. Molecular variance analysis 


(АМОУА) of geographic populations demonstrated that the percentage of variation among the population 


groups of Northeast China, Northwest China, Huang-Huai region, South China, and Southwest China was 


higher than those within each population group. [ Conclusion ] The results revealed that there is less gene 


flow among different geographic populations of L. glycinivorella, which have high-degree whole genetic 


diversity and genetic differentiation. And the observed genetic differentiation primarily occurs among 


population groups of different geographic regions. The geographic distance has important influence on the 


genetic distance of different geographic populations of L. glycinivorella. 


Key words; Leguminivora glycinivorella; mitochondrial СОП gene; geographic population; genetic 


differentiation; gene flow 





大 豆 食 心 虫 Leguminivora glycinivorella (Matsumura), 
别名 豆 葬 过 ,又 名 大 豆 星 莱 蛾 .小 红 虫 . 豆 鞠 虫 等 ,在 
ЖБ ЕЖЕН (Lepidoptera ) 小 卷 蛾 科 
( Olethreutidae) ,是 我 国 大 豆 重要 星 葬 害 虫 ( 史 树 森 
等 , 2014) 。 该 虫 每 年 发 生 1 代 , 并 以 老 熟 幼虫 在 土 
RFRA ,其 食性 单一 , 仅 危害 大 豆 `. 野 生 大 豆 



































和 分 化 等 较 好 的 分 子 标记 ( 李 苹 等 , 2010; 梁 日 霞 
等 , 2011; 李 冉 等 , 2015)。 采 用 CO 本 片段 探讨 大 
豆 食 心 虫 不 同 地 理 种 群 间 的 遗传 分 化 情况 ,可 有 效 
地 研究 大 豆 食心虫 种 内 地 理 分 布 格局 。 本 研究 选用 
线粒体 CO 工作 为 分 子 标记 ,通过 对 中 国 不 同 区 域 大 
豆 食心虫 19 个 地 理 种 群 的 种 群 遗传 结构 和 变异 进 
































和 苦 参 ( 徐 庆 丰 等 , 1965 ) 。 其 幼虫 蚌 入 豆荚 取 食 籽 
粒 , 虫 食 粒 率 轻 者 为 10% ~ 20% , 重 者 可 达 50% 以 
上 ,可 使 大 豆 减 产 20% ~30% ,严重 影响 大 豆 的 产 
量 和 品质 , 且 该 虫 在 我 国 东北 华北 西北 西南 和 华 
东 等 地 均 有 分 布 ( 庞 春 杰 等 , 2012) 。 因 此 ,有 必要 
从 遗传 本 质 上 了 解 全 国 大 豆 食 心 虫 不 同 地 理 种 群 的 
遗传 多 样 性 ,明确 该 害虫 对 不 同 生 态 环境 的 适应 机 
制 ,为 该 害虫 的 综合 治理 提供 理论 依据 。 

国内 外 学 者 对 于 大 豆 食 心 虫 生物 学 特性 ( 徐 庆 
丰 等 ，1965; Sakagami et al., 1985; 史 树 森 等 ， 
2014) ЖЕ МІЗ, 2012) HAR HRE 
等 , 2012; Hu et al., 2013) 、 防 治 技 术 ( 宋 鹏 翔 等 ， 
2014; 032, 2015; Song ей al., 2015) 及 285 
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行 研 究 ,分 析 其 不 同 地 理 种 群 的 遗传 分 化 程度 .遗传 
多 样 性 及 基因 交流 水 平 ,以 期 为 全 面 探讨 中 国 大 豆 
食心虫 不 同 地 理 种 群 间 的 系统 发 育 及 进化 关系 提供 
更 全 面 完整 的 遗传 学 证 据 , 进而 为 该 害虫 的 区 域 性 
发 生 规律 研究 和 防治 策略 制定 等 提供 分 子 生 物 学 基 
础 理论 依据 。 























1 材料 与 方法 


1.1 供 试 材料 

1.1.1 供 试 虫 源 : 本 研究 以 中 国 10 省 (区 ) 共 19 
个 地 理 种 群 的 大 豆 食心虫 老 熟 幼虫 作为 供 试 虫 源 。 
寄主 为 当地 栽培 大 豆 。 将 成 熟 期 的 大 豆 收 割 后 置 于 





























rRNA 基因 全 序列 的 克隆 及 表达 分 析 ( 李冬梅 等 ， 
2012) 等 方面 进行 了 研究 ,但 迄今 有 关 大 豆 食 心 虫 
地 理 种群 的 遗传 结构 及 分 化 方面 的 报道 较 少 ,仅见 王 
红 等 (2014) 采 用 COI 基 因 和 Wang 等 (2015 ) 采 用 COI 
和 Cytb 基因 对 中 国 东 北 地 区 大 豆 食心虫 地 理 种 群 进 
行 了 遗传 多 样 性 及 基因 流 分 析 , 人 研究 结果 显示 大 豆 食 
心 虫 种 群 内 遗传 多 样 性 水 平 处 于 中 低 等 水 平 ,种 群 间 
的 基因 交流 并 未 受到 地 理 距 离 的 影响 。 但 该 研究 仪 
针对 东北 春 大 豆 区 不 同 大 豆 食心虫 种 群 。 

大 豆 食心虫 在 我 国 各 大 豆 产 区 均 有 分 布 ,取样 
范围 和 区 域 代表 性 对 研究 其 地 理 种 群 分 化 至 关 重 
要 。 线粒体 ОМА ( mitochondrial DNA, mtDNA ) 中 的 
CO 了 片段 已 成 为 研究 昆虫 系统 发 育种 群 遗传 变异 
























































铺设 有 塑料 布 的 田间 地 面 或 晒 场 水 泥 地 面 上 ,于 当 
日 中 午 采 集 自主 脱 美 的 老 熟 幼虫, 将 其 浸泡 于 无 水 
乙醇 中 ,并 标注 样品 信息 。 浸 泡 的 第 一 个 月 需 更 换 
无 水 乙醇 3 ~4 次 ,后 置 于 -70Y 超 低温 冰箱 中 保存 
备用 。 各 地 理 种 群 的 具体 采集 信息 及 地 理 群 体 缩写 
代码 详 见 表 1。 
1.1.2 供 试 试剂 : RNase A 酶 .和 蛋白酶 K、DL2000 
DNA Marker 6 x Loading buffer .50 х ТАЕ ЛЕ №, 
Taq PCR Master Міх 扩 增 试剂 盒 和 Ezup 柱 式 动物 基 
因 组 DNA 抽 提 试剂 盒 均 购 自生 工 生 物 工 程 上 海 
(股份 ) 有 限 公司 。 
1.2 大 豆 食心虫 基因 组 ОМА 的 提取 

19 个 大 豆 食心虫 地 理 种 群 的 单 头 老 熟 幼虫 经 
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Жі 中 国 大 豆 食 心 虫 不 同 地 理 种 群 样本 信息 


Table 1 Sampling data of different geographic populations of Leguminivora glycinivorella in China 












































序号 种 群 代码 样本 数 采集 地 点 地 理 坐标 采集 时 间 
No. Population code Number of specimens Collecting locality Geo-coordinates Collecting date 
1 DA 11 广西 都 安 Ои’ ап, Guangxi 23°94'N, 108°09'Е 2013.10 
2 СҮ 12 贵州 贵阳 Guiyang, Guizhou 26°50'N, 106°65'Е 2014.8 
3 DF 8 贵州 大 方 Dafang，Guizhou 27°08'N, 105°35'Е 2015.8 
4 $7, 11 安徽 宿州 Suzhou, Anhui 33°63'№, 116°97'Е 2015.9 
5 XZ 14 江苏 徐州 Xuzhou，Jiangsu 349127У, 117°16'E 2014.10 
6 ҮА 10 陕西 延安 Уап’ ап, Shaanxi 36%60"У, 109°47'Е 2015.10 
7 CZ 10 河北 沧州 Cangzhou，Hebei 38°33'№, 116°83'Е 2013.10 
8 PLD 11 辽宁 普兰 店 Pulandian ，Liaoning 39923"Х, 121°58'Е 2014.9 
9 SY 12 辽宁 沈阳 Shenyang，Liaoning 41948"У, 123°25'Е 2014.10 
10 BX 11 辽宁 本 溪 Benxi，Liaoning 419177Х, 123°45'Е 2015.10 
11 TL 12 辽宁 铁岭 Tieling，Liaoning 42°17'N, 123°50'E 2014.10 
12 СЕ 10 内 蒙古 赤峰 Chifeng, Inner Mongolia 42°28'N, 118°87'Е 2014.10 
13 HL 12 吉林 和 龙 Helong, Jilin 42°32'N, 129°0'E 2014. 10 
14 CC 12 吉林 长 春 Changchun Jilin 43°05'N, 124°18'Е 2015.10 
15 АТ 11 吉林 安 图 Antu, Jilin 43°06'N, 128°53'Е 2014.10 
16 5Н 12 黑龙 江 绥化 Suihua, Heilongjiang 469397Х, 126°57'Е 2014.10 
17 JMS 9 黑龙 江 佳木斯 jiamusi，Heilongjiang 469477Х, 130°18'Е 2014.10 
18 ЛЇТ 10 内 蒙古 扎 兰 屯 Zhalantun, Inner Mongolia 48°0'3, 122°43'Е 2015.10 
19 HH 10 黑龙 江 黑 河 Heihe，Heilongjiang 50°15'№, 127°28'Е 2014.9 


双 蒸 水 漂洗 、 吸 水 纸 吸 干 , 液 氮 研 磨 成 粉末 后 采用 
Елар 柱 式 动物 基因 组 DNA 抽 提 试剂 盒 提取 虫 体 基 
因 组 DNA。 提 取 的 DNA 样品 经 0. 8% 的 琼脂 糖 凝 
胶 电 泳 检测 ,电泳 条 带 需 清晰 明亮 ,合格 样品 保存 于 
-20 冰箱 中 备用 。 

13 大豆 食心虫 线粒体 DNA CO IHA H f PCR 
扩 增 和 产物 测序 

19 个 大 豆 食心虫 地 理 种 群 线粒体 DNA СОП Ж 
因 片 段 PCR 扩 增 上 游 、 下 游 引 物 参 照 Simon 等 
(1994) 的 通用 引物 :TL2-J-3037: 5'-ATGGCAGATT 
АТАТСТААТСС-3”; ТК-№-3785; 5'-СТТТААСАСА 
ССАСТАСТТС-3', 

PCR 扩 增 采用 Тад РСК Master Міх 扩 增 试剂 盒 ， 
反应 体系 为 50 uL: Taq PCR Master Mix 25 uL, 
ddH,O 18 pL,19 个 地 理 种 群 大 豆 食心虫 基因 组 DNA 
模板 (10 pg/mL)2 pL, 上 下 游 引物 (10 mmolL) 各 2 
pL，Mg “(25 mmol/L) 1 pL。 扩 增 反应 条 件 :94C 预 
变性 2 тіп; 94% 变 性 30 в, 52% 退 火 30 s, 72°С 延伸 
1 min， 共 进行 35 个 循环 ;最 后 72°С 延伸 10 шіп, 
PCR 产物 5 pL 于 0. 8% 的 琼脂 糖 凝 胶 进行 电泳 检 
测 , 凝 胶 成 像 分 析 系 统 观 察 并 拍照 ,确定 扩 增 成 功 的 
委托 华 大 基因 科技 有 限 公司 测序 。 

1.4 数据 分 析 
将 测序 所 得 的 19 个 大 豆 食 心 虫 地 理 种 群 线 粒 












































Ж DNA CO 基因 序列 采用 Chromas 软件 读 取 并 观 
察 峰 值 , 通 过 MEGA6. 0 (Tamura et al., 2013) 软件 
中 ClustalW 进行 多 重 序列 比 对 ,去 除 序列 两 端 测 序 
误差 较 大 的 区 域 ,将 处 理 后 的 序列 在 NCBI 网 站 上 
进行 BLAST 比 对 、 确 认 。 应 用 MEGA6.0 软件 分 析 
碱 基 组 成 .保守 位 点 、 变 异 位 点 、 简 约 信息 位 点 、 自 疹 
位 点 、 转 换 和 站 换 碱 基 对 数 以 及 计算 单 倍 型 间 的 遗 
传 距 离 ; УЛ ЗЕҢ Суда pomonella (L. ) К) СОП 
基因 片段 (GenBank 登录 号 : JX407107. 2 ) 为 外 群 ， 
采用 最 大 似 然 法 (maximum likelihood, М) 构建 不 
同 地 理 种 群 的 大 豆 食心虫 单 倍 型 聚 类 树 ,并 对 各 分 
文 的 置信 和 度 自 展 法 (Bootstrap ) 进行 1 000 次 重复 检 
验 ; 采 用 不 加 权 配 对 算数 法 (unweighted pair-group 
method with arithmetic means, ОРСМА ) 以 不 同 地 理 
种 群 间 遗 传 距离 矩阵 构建 不 同 地 理 种 群 间 系 统 发 育 
树 。 利 用 DnaSP 5.0 软件 (Librado апа Rozas, 2009) 
分 析 19 个 地 理 种 群 内 CO 基因 序列 单 倍 型 数目 和 出 
现 的 频率 、 单 倍 型 多 样 度 Hd ЕЕ НЕЕ Pi It 
酸 平均 差异 数 K、 基 因 流 Nm、 遗 传 分 化 系数 Gst、 固 定 
系数 及 上 种 群 间 核 昔 酸 平均 差异 数 Кху IAE REX 
Ж Dxy 等 分 子 遗 传 学 系数 并 进行 中 性 检验 (Tajima， 
1989) 。 依 据 经 纬度 计算 各 采集 地 之 间 的 地 理 距离 ， 
使 用 GenAlEx6. 41 软件 中 Mantel 检验 (Peakall апа 
Smouse , 2006 ) 检验 不 同 地 理 种 群 间 遗传 距离 与 地 理 
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距离 自然 对 数值 的 相关 性 。 单 倍 型 网 络 图 由 
Network5.0 软件 (Polzin апа Daneschm , 2003) 绘制 。 

19 个 地 理 种 群 按照 其 所 在 区 域 共 分 为 5 个 地 理 
种 群 组 :都 安 (DA ) 属 华南 组 ;延安 (YA) 属 西北 组 ; 贵 
阳 (GY) 和 大 方 (DF) 属 西南 组 ;宿州 (SZ) RIM (XZ) 
及 沧州 (CZ) 属 黄 淮 组 ;普兰 店 (PLD) 沈阳 (SY) Ж 
溪 (BX) ШЕТІ), (СЕ) 、 扎 兰 屯 (ZLT) 绥化 
(SH) JER Hr (JMS). (HH) MJ (HL), KE 
(СС) 和 安 图 (AT) 划 分 为 东北 组 。 根 据 以 上 分 组 ,应 
用 Arlequin3.5.1.2 软件 (Excoffier et al., 2005 ) 进行 
АМОУА 分 子 变 异 分 析 以 评估 大 豆 食心虫 的 遗传 分 
化 情况 ,同时 应 用 DnaSP 5.0 软件 分 析 5 组 大 豆 食 心 
虫 地 理 种 群 的 基因 流 Nm 和 固定 系数 Fst。 




















2 结果 


2.1 碱 基 组 成 及 序列 变异 

CO 本 基因 片段 预期 扩 增 长 度 为 750 bp ,经 多 重 
序列 比 对 剪 切 后 ,本 研究 最 终 获得 19 个 大 豆 食心虫 
地 理 种 群 线粒体 CO 了 基因 长 度 为 723 bp 的 序列 片 
段 。 本 研究 共 获 得 208 条 序列 , 均 没有 碱 基 的 缺失 
和 插入 ,保守 位 点 662 个 (91. 56% ) ;变异 位 点 61 个 
(8.44% ) ,其 中 简约 信息 位 点 39 1-(63.93%),Н 
裔 位 点 22 个 (36. 07% ) 。 核 苷 酸 蔡 换 分 析 表 明 , 转 
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图 1 大 豆 食 心 虫 不 同 地 理 种 群 线粒体 CO 开 基 因 各 单 倍 型 的 中 介 


Haplotype network of different geographic populations of Leguminivora glycinivorella а оп mitochondrial СОП gene 


换 的 位 点 数 55 个 ,A 与 С 之 间 的 转换 为 16 个 ,C 与 
T 之 间 的 转换 为 39 T; RERA 6 个 ; 另 
有 1 个 C/G/T 和 1 个 А/С/Т 三 碱 基 变异 位 点 ,总 
体 序 列 变异 中 转换 多 于 颠 换 ,符合 亲缘 关系 较 近 的 
2. Ж 
1994) 。 所 测序 列 CO 基因 片段 组 
成 中 ,A + 了 平均 含量 为 (75. 60% ) ,明显 高 于 C +C 
含量 (24. 40% ) ,表现 出 明显 的 碱 基 偏 向 性 ,符合 昆 
虫 线粒体 基因 序列 碱 基 组 成 特点 (Liu and 
Вескепраһ, ，1992 ) 。 
22 遗传 多 样 性 和 单 倍 型 分 析 

在 19 个 大 豆 食心虫 地 理 种 群 中 共 获 得 57 种 单 
倍 型 ,分别 记录 为 H1 - H57(KX119989 - KX120045 ) 。 
应 用 Kimura-2-Parameter 模型 计算 单 倍 型 之 间 的 遗传 
距离 ,结果 显示 不 同 单 倍 型 遗传 距离 在 0.0010 ~ 
0. 0420 之 间 , 平 均 遗 传 距离 为 0.0090。 单 倍 型 H3 出 
现 频 率 最 高 (70/208) , 占 所 有 检测 个 体 的 33.65% ,出 
现在 所 检测 的 18 个 地 理 种 群 中 ;其 次 出 现 频率 较 高 
的 单 倍 型 是 НА (25/208), 占 所 有 检测 个 体 的 
12. 01% ,出 现在 所 检测 的 11 个 地 理 种 群 中 。 单 倍 型 
网 络 中 介 图 显示 (图 1) , 单 倍 型 H3 位 于 单 倍 型 网 络 
中 介 图 的 中 心 ,其 他 单 倍 型 围绕 其 分 布 , 单 倍 型 H19, 
H20, H21, H27 和 H28 与 其 他 单 倍 型 相对 遗传 距离 
较 远 , 且 与 其 他 单 倍 型 有 明显 分 支 。 
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НІ - H57: 在 19 个 大 豆 食 心 虫 地 理 种 群 中 定义 的 单 倍 型 Haplotypes of 19 geographic populations of L. glycinivorella. 下 同 Тһе same below. 圆 面积 











代表 单 倍 型 出 现 频率 。Circle areas аге proportional to haplotype frequencies. 
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根据 DnaSP5.0 软件 分 析 结 果 表 明 ( 表 2) ,大 豆 
食心虫 19 个 地 理 种 群 间 存 在 一 定 的 遗传 差异 。 总 
群体 单 倍 型 多 样 度 Hd =0. 8522 , 6 0 ЖЕТЕРІ = 
0. 0065 , 核 苔 酸 平均 差异 数 K=4.6995。 各 地 理 种 群 
单 倍 型 数量 范围 在 2 ~ 10 个 之 间 , 长 春 ( CC) 种 群 单 
倍 型 数量 最 丰富 ,在 检验 个 体 中 定义 了 10 种 单 倍 
型 。 各 地 理 种 群 的 单 倍 型 多 样 度 Hd 范围 为 








0. 1818 ~ 0. 9546 ,其 中 宿州 (SZ) 种 群 单 倍 型 多 样 度 
最 低 , 长 春 (CC) 种 群 单 倍 型 多 样 度 最 高 。 从 
Tajima’s D 检测 结果 可 以 得 出 ,19 个 地 理 种 群 的 中 
性 检测 结果 除 徐州 (XZ) 外 均 不 显著 , 上 且 总 群体 的 
Tajima’ s D 值 为 负 值 ,结果 不 显著 (0.10 >P > 
0.05) ,表明 19 个 地 理 种 群 的 大 豆 食心虫 总 体 上 在 
较 近 的 历史 时 期 未 经 历 明 显 的 种 群 扩张 。 

































































R2 大 豆 食心虫 不 同 地 理 种 群 CO 基因 单 倍 型 多 样 度 , 核 昔 酸 多 样 性 分 析 及 Tajima 's D 中 性 检验 
Table 2 СОП haplotype diversity, nucleotide diversity and Tajima’ s D test of 
different geographic populations of Leguminivora glycinivorella 
SA нл 单 倍 型 多 样 度 核 车 酸 多 样 度 ” 核 车 酸 平均 差异 类 ыы 
种 群 代码 РАВ ои о оО 
i Haplotype Haplotype Nucleotide Average number 中 性 检测 ЕС 
Population Г А 4 ; | ў К И ШЕ Statistical 
distribution diversity diversity of nucleotide differences Tajima’ s D Км 
code significance 
һ На Pi K 
DA Н16(1) Н19(2) Н20(5)Н21(3) 0. 6909 0. 0065 4.6909 -0.7498 Р>0.10 
GY H3 (1 ) H23 (2) Н24(5)Н27(3)Н28(1) 0.7879 0.0170 12.273 1.0183 P >0.10 
52 Н3(10) Н43(1) 0.1818 0. 0003 0.1818 -1. 5334 P >0.10 
H3(2)H14(1)H16(3) 
YA 0.888 0. 0035 2.5111 – 1.1028 Р>0.10 
Н35(2) Н40(1)Н52(1) 9 5 
H3(5)H4(1)H14(1) 
CZ 0.7778 0.0023 1.6667 -1.3827 P >0.10 
H16(1)H17(1)H18(1) > 
н3(4)Н4(2)Н15(1) 
СЕ 0.8667 0.0002 1.5556 -1.1028 P >0.10 
Н16(1) Н25(1)Н26(1) 42 
PLD Н3(7)Н4(1)Н14(1)Н40(1)Н41(1) 0.6182 0.0015 1.0546 - 1.4646 P >0.10 
СС H3 (3) HA CT) HS CT) H8 (1) H? НЮ 0.9546 0. 0041 2.9697 -1.5334 P >0.10 
H11(1)H2(1)H13(1)H14(1) 
HH Нано 0.9111 0.0023 1.6667 -1.3827 P >0.10 
H31(1)H32(1)H33(1) 
AT H12) H2 C1) H3 (2) HACS) 0.9091 0.0031 2.2364 -1.5916 0.10>Р>0.05 
Н5(1)Н6(1)Н7(1) 

DF H3 (3) H22 (1)H23(1)H24(3) 0.7857 0.0049 3.5714 -0.8552 Р>0.10 
JMS HS ШО ТУНЫ 30 0.9444 0.0033 2.3889 1.4646 P>0.10 
Н37(1)Н38(2)Н39(1) ; i i i 24 
ВХ НЗ(1)НІ4(4) Н38(6) 0.6182 0.0024 1.7455 — 1.3298 Р>0.10 
SY H3(6)H4(1)H14(2)H42(1)H43(1)H44(1) 0.7576 0.0023 1.6818 -1.0740 Р>0.10 
5Н H3(2)H4(6)H45(1)H46(1)H47(1)H48(1) 0.7576 0.0021 1.5303 -1.1285 Р>0.10 

H3(4)H4(1)H9(1)H16(1)H25(1) 

TL 0.893 0.0024 1.7121 —1.41244 P>0.10 
Н35(1)Н40(2)Н49(1) <р” z 

XZ НЗ (12) Н50(1) Н51(1) 0. 2747 0. 0008 0.5714 -1.79759 Р <0.05 

ZLT HI CLHA JHS еи 0. 9333 0. 0039 2.8444 -0.48329 P >0.10 
Н55(1) Н56(1) Н57(1) 

HL та (2) HECH) H1S (19H21) 0. 8939 0. 0028 2.0303 -0.86572 P >0.10 
Н30(1) Н34(1) НЗ5(1) 

总 种 群 
Total 0.8522 0.0065 4.6995 -1.74916 0.10>Р>0.05 
populations 
种 群 代码 信息 见 表 1。Population code information see Table 1. 下 同 Тһе same below. 
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2.3 系统 发 育 分 析 

CO 了 基因 单 倍 型 ML 树 显示 (图 2), 单 倍 型 
H19, H20 #1 Н21 是 都 安 (DA ) 特有 单 倍 型 聚 为 一 
支 , 单 倍 型 H27 和 H28 为 贵阳 (GY) 特有 单 倍 型 亦 
限 为 一 支 , 贵 阳 和 都 安 的 特有 单 倍 型 与 其 他 地 理 种 
群 间 和 遗传 距离 相差 较 大 。 在 所 构建 单 倍 型 系统 发 育 
树 中 ,尽管 有 距离 较 近 地 理 种 群 的 单 倍 型 优先 聚 在 
一 个 分 文 [ 如 黑河 (HH) 的 H33 和 绥化 (SH) 的 
H48 ] ,但 也 有 地 理 距 离 相距 较 远 的 单 倍 型 聚 在 一 起 
的 情况 [如 扎 兰 屯 (ZLT) 的 H56 和 安 图 (AT) 的 
H6 ] ,是 各 分 支 间 置 信 度 较 低 , 说 明 这 些 单 倍 型 之 间 
差异 较 小 ,不 能 形成 可 靠 的 分 支 。 基 于 大 豆 食 心 虫 
不 同 地 理 种 群 间 的 遗传 距离 构建 的 UPGMA RAH 
显示 (图 3) ,其 分 支 情况 与 单 倍 型 系统 发 育 树 情况 
一 致 。 即 贵阳 (GY) 和 都 安 ( DA) 种 群 单独 形成 分 
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文 ,形成 具有 较 大 遗传 差异 性 地 理 种 群 ,与 其 他 地 理 
种 群 分 化 明显 。 
2.4 遗传 分 化 分 析 

大 豆 食 心 虫 各 地 理 种 群 总 遗传 分 化 系数 Gst = 
0. 1186 ,总 国定 系数 Fst = 0.4619, 总 基因 流 Nm = 
0. 29 ;大 豆 食 心 虫 东北 组 、 华 南 组 西南 组 西北 组 及 
黄 淮 组 5 组 间 基 因 流 Nm =0.13 ,固定 系数 Fst = 
0.6498。 各 地 理 种 群 间 核 昔 酸 差异 数 Kxy 在 
0. 6364 ~ 23. 8485 È |н], М F R Iy XE Dxy 在 
0.0016 ~ 0.0290 之 间 , 遗传 分 化 系数 Ся 在 
-0.2134 ~ 0. 3924 之 间 , 固 定 系数 Fst 在 -0.0547 ~ 
0. 8744 2 |н], ABZ (DA) 和 贯 阳 (GY) 种 群 与 其 他 
地 理 种 群 间 的 Fst 值 在 0. 1606 ~ 0. 8744 之 间 , 遗 传 
分 化 水 平 较 高 ,其 余 17 个 地 理 种 群 间 遗 传 分 化 程度 
则 相对 较 低 ( 表 3) 。 
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最 大 似 然 法 (ML) 法 构建 大 豆 食心虫 不 同 地 理 种 群 线粒体 CO I ЖЕГЕ EE Ж ОУ 








Fig. 2 Maximum likelihood (ML) phylogenetic tree of mitochondrial СОП gene haplotypes from 


different geographic populations of Leguminivora glycinivorella 
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ЗИ СОП ЗЕ УУК. СОП gene sequence of Cydia pomonella was used as the outgroup. 
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ХАН (ЖЕН 西南 组 西北 组 及 黄 淮 组 进 
行 AMOVA 分 析 后 显示 ,组 间 变 异 占 方差 比率 的 
57. 01% ,组 内 种 群 间 变 异 占 方差 比率 的 2 58% ,种 
群 内 变异 占 方差 比率 的 40. 41% ,表明 遗传 变异 主 
要 来 自 组 间 ( 表 4)。 采 用 Kimura-2-Paramter 模型 分 
析 大 豆 食 心 虫 不 同 地 理 种 群 遗传 距离 ,发 现 不 同 地 
理 种 群 遗 传 距离 在 0.0005 ~ 0. 0567 之 间 , 东北 组 
12 个 地 理 种 群 间 的 遗传 距离 在 0. 0018 ~ 0. 0038 之 
间 ( 表 3)。 除 都 安 (DA) 和 贵阳 (GY) 与 其 他 地 理 种 
群 之 间 的 遗传 距离 较 大 外 ,其 余 均 较 低 。 东 北 组 各 
地 理 种 群 间 遗 传 距离 与 地 理 距 离 的 Mantel 相关 性 
分 析 表 明 ( 图 4: B) ,各 地 理 种 群 间 基因 交流 受 地 理 
距离 的 影响 较 小 (r = -0.02, P=0.56 >0.05); 但 











在 全 国 范围 内 ,19 个 地 理 种 群 间 遗传 距离 与 地 理 距 
离 的 Mantel 相关 性 分 析 表明 (图 4: A) ,大 豆 食 心 忠 
各 地 理 种 群 间 地 理 距 离 与 遗传 分 化 程度 存在 相关 性 
(r =0. 394, 已 =0.01 <0. 05) ,说 明 在 较 大 区 域 范 围 
内 ,地理 距 离 是 影响 大 豆 食 心 虫 不 同 地 理 种 群 间 遗 
传 距离 和 基因 交流 的 重要 因素 之 一 。 




















3 讨论 


祖 元 刚 等 (1999 ) 认为 一 个 物种 群体 内 部 的 遗 
传 多 样 性 能 反映 此 物种 对 环境 变化 的 适应 能 力 。 本 
试验 研究 结果 显示 ,19 个 大 豆 食心虫 地 理 种 群 总 群 
体 的 单 倍 型 多 样 度 友 =0. 8522 , 共 57 种 单 倍 型 , 表 
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图 3 UPGMA 法 构建 大 豆 食心虫 19 个 地 理 种 群 线粒体 СО IERRA 





Fig. 3 UPGMA phylogenetic tree of mitochondrial СОП репе haplotypes from 19 


geographic populations of Leguminivora glycinivorella 
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图 4 基于 线粒体 COI 基 因 大 豆 食心虫 不 同 地 理 种 群 间 遗 传 距离 与 地 理 距 离 对 数 的 Mantel 检验 


Fig. 4 Mantel test between genetic distance апа the logarithm of geographic distance among different 


geographic populations of Leguminivora glycinivorella based on mitochondrial СО 1 gene 


A: 19 个 地 理 种 群 Nineteen geographic populations; B: 东北 组 Northeast China population group. 
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表 4 ”大豆 食 心 虫 不 同 地 理 种 群 组 间 遗 传 变异 的 分 子 变异 分 析 
Table 4 AMOVA analysis of different geographic population groups of Leguminivora glycinivorella 
变异 来 源 自由 度 平方 和 方差 组 分 变异 百分率 固定 指数 

Source of variation df Sum of squares Variance components Percentage of variation Fixation index 
组 间 Among groups 4 213.86 1.78 Иа 57.01 Есі -0.57 (Р<0.01) 
组 内 种 群 间 Among populations within а group 14 30.02 0.08 Vb 2.58 Езс -0.06 (Р<0.05) 
种 群 内 Within а population 189 238.43 1.26 Үс 40.41 Fst =0.60 (Р<0.00001) 


Уа, Vb, Үс: 方差 组 分 的 数量 Number of variance components. 
明 大 豆 食心虫 种 群 具有 较 高 的 单 倍 型 多 态 性 ,整体 


的 遗传 多 样 性 相对 较 高 ; 杨 哲 等 (2015 ) 采用 СО I 
基因 研究 亚洲 玉米 蜡 Ostrinia furnacalis 种 群 间 遗传 


其 他 种 群 的 Fst 值 在 0. 1185 ~ 0.3037 之 间 , 大 多 处 
于 中 等 遗传 分 化 水 平 , 且 大 方 (DF ) 种 群 在 两 个 系统 
树 中 均 与 贵阳 (CY ) 种 群 不 在 同一 分 文 上 , 究 其 原 











多 样 性 , 共 检 测 出 48 种 单 倍 型 , 认为 亚洲 玉米 旦 对 
不 同 寄主 植物 环境 具有 较 强 的 适应 能 力 和 较 高 的 遗 
传 变 异 潜力 ,与 本 试验 结果 类 似 ,同时 在 本 研究 中 发 
现 19 个 地 理 种 群 的 单 倍 型 多 样 度 与 核 音 酸 多 样 度 
均 存 在 一 定 差异 ,体现 出 不 同 的 遗传 多 样 性 水 平 。 
孙 党 等 (2013) 研 究 认 为 ,可 能 是 生态 栖息 环境 的 多 
ЖЕЕ ‚ж ВИД Ду ЕШ Oedaleus infernalis PREEN E 
富 的 遗传 多 态 性 ; 李 春 选 等 (2003 ) 研究 表明 ,中 华 
K ЭЙ Ж Ж ОЕШ Shirakiacris 
shirakii, 18 Haplotropis brunneriana №. 55 Zn 1А, Їй 
Atractomorpha sinensis 因 其 对 栖息 环境 的 依赖 程度 ， 
体现 出 不 同 的 遗传 多 样 性 水 平 ,由 此 推测 大 豆 食 心 
虫 不 同 地 理 种 群 遗传 多 样 性 水 平 不 同 的 原因 可 能 是 
因为 其 成 虫 飞行 能 力 较 弱 , 地 理 屏 障 如 峡谷 、 山 脉 等 
会 有 效 限 制 大 豆 食心虫 的 扩散 和 迁移 ,同时 野生 大 
豆 分 布 及 栽培 大 豆 区 域 范围 等 寄主 条 件 也 是 其 发 生 
分 布 的 限制 因素 。 

在 19 个 大 豆 食心虫 地 理 种 群 获得 的 57 种 单 倍 
型 中 , 单 倍 型 H3 НА, H14 和 H16 至 少 被 6 个 种 群 
共享 ,为 主体 单 倍 型 ; 单 倍 型 H19, H20, H21 和 
H27 等 单 倍 型 为 某 个 地 理 种 群 所 特有 ,为 独 享 单 倍 
型 。 独 享 单 倍 型 和 共享 单 倍 型 的 出 现 则 说 明 各 地 理 
种 群 在 存在 着 一 定 基因 交流 的 同时 ,也 具有 一 定 程 
度 上 的 遗传 分 化 ( 孙 昌 等 ,2013; ІРІҢ, 2016). 
大 豆 食 心 虫 总 群体 Tajima’ s D AEN -1.74916 Н. 
检验 结果 不 显著 ,表明 序列 在 进化 上 符合 中 性 模型， 
大 豆 食心虫 总 体 在 较 近 的 历史 时 期 内 未 经 历 群 体 扩 
张 ,种 群 大 小 保持 相对 稳定 的 状态 (Harpending ей 
al., 1998; 王 兴 亚 和 许 国庆 , 2014); 

Fst 是 衡量 遗传 分 化 程度 的 一 个 重要 参数 ,西南 
的 贵阳 (GY) 和 华南 组 的 都 安 ( DA ) 种 群 与 东北 
旦 西北 组 及 黄 淮 组 间 的 Fst 值 达到 0.2441 以 上 ， 
遗传 分 化 程度 较 大 。 但 西南 组 的 大 方 (DF ) 种 群 与 
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因 可 能 是 因为 大 方 (DF) 种 群 样本 的 采样 地 位 于 贵 
州 省 西北 部 ,地 处 云贵 高 原 向 黔 中 原 丘 陵 过 渡 的 斜 
坡地 带 , 属 北 亚热带 湿润 季风 气候 ,相对 为 高 纬度 
(27°08'\) 高 海拔 (1 512 m) 地 区 ;而 贵阳 种 群 采 集 
地 位 于 贵州 省 中 部 ,地 处 丘 原 盆地 地 区 , 属 亚热带 湿 
润 温 和 型 气候 ,其 纬度 为 26"50'N ,海拔 为 1 134. 2 
m ,纬度 与 海拔 均 低 于 大 方 县 。 相 关 研 究 表明 ,海拔 
高 度 每 增高 100 mm ,气温 降低 0.5 ~0. 6Y ,光照 增加 
4% ~5% , 即 大 方 地 区 的 光照 时 间 长 于 贵阳 地 区 , 温 
度 则 低 于 贵阳 地 区 ,其 温度 环境 更 趋 近 于 高 纬度 地 
区 。 史 树 森 等 (2014 ) 人 研究 表明 温度 对 大 豆 食心虫 
卵 的 发 育 速率 、 幼 虫 存 活 率 、 幼 虫 期 营养 积累 均 有 影 
响 , 且 大 豆 食心虫 属 “ 长 日 照 昆 虫 ” ,因此 笔者 认为 
贵阳 (GY) 与 大 方 (DF) 种 群 间 产 生 一 定 遗 传 分 化 的 
原因 可 能 是 二 者 所 处 气候 、 生 态 环境 存在 较 大 差异 
所 致 。 

王 红 等 (2014) 和 Wang 等 (2015 ) 对 中 国 东北 地 
区 大 豆 食 心 虫 地 理 种 群 进行 了 遗传 多 样 性 研究 , 结 
果 表 明 大 豆 食心虫 种 群 遗 传 分 化 程度 较 小 , 基因 交 
流 未 受 地 理 距 离 影 响 。 而 笔者 研究 了 全 国 各 大 豆 产 
区 的 19 个 大 豆 食 心 虫 种 群 ,结果 表明 ,总 群体 固定 
系数 Fst =0. 4619 ,说 明 大 豆 食 心 虫 不 同 地 理 种 群 间 
存在 很 大 程度 的 遗传 分 化 ; 总 群体 基因 流 Nm = 
0. 29 ,表明 各 地 理 种 群 间 基 因 流 水 平 很 低 , 可 能 由 于 
遗传 漂 变 发 生 了 分 化 ( 张 利 娟 等 , 2012) ІНІМ, Ж 
ДЕН 西北 组 西南 组 、 黄 淮 组 、 华 南 组 5 组 间 数 据 分 
析 结 果 显 示 固 定 系数 Fst = 0. 6498 , 基因 流 Nm = 
0. 13 ,表明 组 间 遗 传 分 化 程度 较 大 。 结 合 AMOVA 
分 析 和 Mantel 检验 结果 可 知 , 大 豆 食 心 虫 遗传 变异 
主要 来 自 组 间 , 且 大 豆 食 心 虫 种 群 间 的 地 理 距 离 与 
遗传 分 化 程度 存在 显著 相关 性 ,基因 交流 受 地 理 距 
离 影 响 。 但 在 本 研究 中 ,如 果 将 农业 生态 环境 相对 
一 致 的 东北 组 12 个 地 理 种 群 单独 分 析 , 其 遗传 距离 
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在 0. 0018 ~ 0. 0038 之 间 ,遗传 距离 相对 较 小 ;东北 组 
地 理 种 群 的 大 多 数 单 倍 型 在 ML 单 倍 型 系统 发 育 树 
中 都 处 于 较为 混杂 的 分 布 格局 ,未 显示 出 明显 的 地 
理 谱系 , 且 东 北 组 的 地 理 种 群 在 UPGMA 聚 类 树 中 
聚 成 一 文 , 其 地 理 种 群 间 Fst 值 为 -0.5470 ~ 
0. 1517 之 间 ,遗传 分 化 程度 较 低 ; 将 东北 组 各 地 理 
种 群 间 遗传 距离 与 地 理 距离 的 Mantel 相关 性 分 析 
显示 ,种 群 间 基 因 交 流 受 地 理 距离 的 影响 较 小 ,结果 
与 王 红 等 (2014) 和 Wang 等 (2015 ) 研究 的 大 豆 食 心 
虫 中 国 东北 10 个 地 理 种 群 遗 传 多 样 性 分 析 结 果 类 
似 。 可 见 , 只 有 当地 理 距 离 足以 产生 明显 环境 异 质 
性 的 时 候 , 才 能 导致 大 豆 食 心 虫 地 理 种 群发 生 遗 传 
分 化 。 

对 植 食 性 昆虫 来 说 ,异域 分 布 的 昆虫 产生 种 下 
分 化 是 由 地 理 和 寄主 的 双重 作用 决定 的 (Pashley， 
1986)。 在 本 研究 中 ,贵阳 和 都 安 所 处 地 理 纬 度 
(23°94' ~26°50'N) 明 显 低 于 其 他 地 区 (33°63'N 以 
北 ) , 属 南方 多 作 大 豆 区 ,大 豆 栽培 品种 和 耕作 制度 
多 样 , 与 其 他 种 群 所 处 地 区 的 地 理 位 置 海拔 气候 
及 大 豆 种 植 品种 、 耕 作 制 度 均 存 在 显著 差异 。 张 健 
等 (2015 ) 研究 表明 ,环境 因素 对 北美 五 倍 子 是 
Melaphis rhois 的 遗传 分 化 影响 较 大 ,地 理 隔 离 能 够 
阻碍 蚜虫 基因 交流 ,各 地 理 种 群 间 分 化 显著 ;Tree 
和 Walter(2009) 人 研究 表明 , 榕 属 Ficus 不 同 品 种 寄主 
对 母 管 萄 马 属 Mesothrips 获 马 的 行为 特性 和 生活 史 
具有 一 定 影响 。 因 此 ,综合 分 析 单 倍 型 在 中 介 网 络 
图 和 实际 地 理 种 群 中 的 分 布 情况 、 单 倍 型 系统 进化 
树 和 以 不 同 地 理 种 群 间 遗传 距离 构建 系统 发 育 树 的 
分 析 结 果 来 看 ,地 理 距 离 .纬度 海拔 的 差异 .栽培 耕 
作 制 度 及 大 豆 种 植 品 种 等 差异 均 可 能 是 造成 大 豆 食 
心 虫 不 同 地 理 种 群 产生 遗传 分 化 的 重要 原因 。 

本 研究 仅 初步 明确 了 中 国 大 豆 食 心 虫 不 同 地 理 
种 群 间 的 遗传 分 化 及 其 与 环境 的 关系 ,但 因 所 选 序 
列 相 对 较 短 ,个别 地 区 样本 数量 较 少 ,所 获得 的 遗传 
学 数据 尚 不 够 全 面 。 因 此 在 下 一 步 的 研究 中 ,应 增 
加 试验 采样 地 及 样本 数量 和 所 选 分 子 标 记 测 序 长 
E ,并 配合 其 他 分 子 标记 方法 ,以 便 深 入 探讨 中 国 大 
豆 食心虫 不 同 地 理 种 群 间 遗 传 分 化 状况 及 其 形成 
机 制 。 
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